Far-field spectral characterization of conical emission and filamentation in Kerr media 
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By use of an imaging spectrometer we map the far-field (9 — A) spectra of 200 fs optical pulses 
that have undergone beam collapse and filamentation in a Kerr medium. By studying the evolution 
of the spectra with increasing input power and using a model based on stationary linear asymptotic 
wave modes, we are able to trace a consistent model of optical beam collapse high-lighting the 
interplay between conical emission, multiple pulse splitting and other effects such as spatial chirp. 

PACS numbers: 190.5940, 320.2250 



I. INTRODUCTION 

Filamentation and, in general, collapse, of high 
power laser pulses in transparent media has attracted 
significant attention ever since its prediction [ESQ 
and observation Q and, due to the complexity of the 
involved phenomena, continues to be a hotly debated 
topic. There are many reasons for this interest ranging 
from application possibilities such as localizing ultra- 
short laser pulses over long propagation distances [3, 
white-light laser sources for parametric amplification |y| 
or spectroscopy 0, H] or the formal analogy with the 
equations that describe nonlinear wave collapse in other 
systems, such as Bose-Einstein condensates [jf , that are 
experimentally less accessible. Numerical investigation 
of optical wave collapse is usually carried out starting 
from the nonlinear Schrodinger equation (NLSE) that 
describes the spatial evolution of a beam considering 
diffraction and a self-focusing term that originates from 
the real part of the third-order medium Kerr nonlinearity 
(nv). Such an equation predicts the formation of an 
unstable 2D stationary solution, the so called Townes 
profile [jj, that, if perturbed, will either diffract or 
undergo catastrophic collapse. However, ultra-short 
laser pulses do not, in general, follow this behavior [To| : 
the collapse is arrested by other effects such as pulse 
lengthening due to normal group-velocity dispersion 
(GVD) or plasma defocusing and an apparently sta- 
tionary propagation regime (filament) is reached. Thus 
the NLSE may be accordingly modified to account for 
space-time coupling and also for plasma generation and 
defocusing. These modified equations have proved to 
be able to describe many of the phenomena associated 
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with pulse filamentation such as conical emission (CE) 
RlJ, super-continuum (SC) generation, pulse steepening 
[l2| and splitting 01 . One of the main features that 
is emerging is the importance of space-time coupling. 
Indeed, there exist regimes characterized simultaneously 
by ultra-tight focusing and ultra-short pulse length in 
which the nonlinearity couples the spatial and temporal 
dynamics following an underlying geometry dictated by 

the modulational instability (MI) gain profile nana, in 

such cases it is preferable to avoid space-time separation 
and refer directly to the geometrical properties of the 
space-time environment. A measure of the importance 
of these new concepts is given, for exam ple , by the 
discovery of nonlinear X- waves 0, 0, 0, 0] , chaotic 
spatio-temporal fragmentation due to space-time MI 
Rjl, red solitons [2(|, X- waves in pulse filamentation 
|2l| and stationary conical waves supported by nonlinear 
losses 22] . The complexity of these issues requires 
a careful examination of the experimental methods 
employed, as space-time coupled phenomena should be 
investigated with adequate instruments. 
We note that the major part of laser physics diagnostics 
is based on the separation of spatial and temporal effects 
leading, for example, to the widely used concepts of 
carrier spatial and temporal frequencies, beam walk-off 
and group velocity, diffraction and dispersion etc., i.e. 
either purely spatial or purely temporal quantities. 
This conceptual space-time division is also reflected in 
the standard experimental characterization methods 
also adopted for the study of beam filamentation. We 
have on the one hand, near and far field imaging that 
give information on the spatial profiles but ignore 
the temporal profile. On the other hand, temporal 
auto-correlation traces or frequency resolved optical 
gating (FROG) give a precise measurement of the 
temporal profile at the medium output of the whole 
(or a particular portion of the) filament [Ijj, |23|, 124) . 
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thus losing any information regarding the (transverse) 
space-dependent temporal profile. The limitations of 
FROG measurements have been p artly overcome with 
the so called SPIDER technique |25j that gives the 
temporal profile across one spatial dimension l2ol and 
by use of a particular 3D optical gating method |27l l28l ] 
that provides the full spatio-temporal intensity profile 
of femtosecond pulses and has indeed allowed a full 
dimensional intensity space-time characterization of 
filaments in water [23 and of X- waves in media [28| . 
Although extremely powerful, this technique requires 
the use of two separate, high power, synchronized laser 
sources, the first providing the pulse under investigation 
and the second a pulse that must have an appreciably 
shorter time duration in order to guarantee a high 
(temporal) resolution. Moreover, this method is based 
on a second-order nonlinear up-conversion process: due 
to the large (spatial and temporal) bandwidths involved 
it is necessary to resort to very thin (20 /im or less) 
conversion crystals which in turn lead to low output 
powers and the necessity to integrate over many pulses. 
Here we propose another possible experimental in- 
vestigation tool, namely far-field spectral (6 — A) 
characterization, that may give readily accessible details 
of the spatio-temporal pulse profile in single shot acquisi- 
tion. This method can be used as a general investigation 
tool for phenomena that involve space-time coupling 
and we show its application to pulse filamentation in 
normally dispersive Kerr media. To the best of our 
knowledge there are only a few papers in literature 
that show the angular (9 — A) spectra of Kerr-induced 
filament s jlOl l30l l31j without actually paying attention to 
the details these may contain and certainly no systematic 
study in this sense has been carried out before. The 
data are tentatively explained with the aid of a simple 
interpretation model, based on a linear description of 
the stationary asymptotic light wave modes that act as 
attractors during the pulse evolution [32| . 



II. SURVEY OF BEAM COLLAPSE AND 
FILAMENTATION OF FS PULSES 

In this work we have used a normally dispersive Kerr 
medium with a fixed total length and the only variable 
parameter was the input pulse energy. Therefore, we 
give here a brief overview of the expected pulse evolution 
for increasing input energies based on the existing 
literature. 

For very low powers no self focusing (SF) occurs and 
the beam behaves linearly. However, with increasing 
power the threshold for catastrophic self-focusing is 
reached and the beam will start to spatially contract. A 
simple expression for the critical power for cw gaussian 
beams is given by |33| 
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FIG. 1: Output beam diameter for varying input pulse en- 
ergy. All measurements reported in this work were performed 
in the energy interval 2-5.6 /iJ where single filamentation was 
observed. 



Pcrit — q (1) 

where the nonlinear refractive index is given as a func- 
tion of the pulse intensity / by n = no + n^I- If the 
medium dispersion is also accounted for then, in gen- 
eral, the threshold for SF will be higher for short pulses 
than for longer ones due to the damping effect of nor- 
mal GVD and eq. Q should be accordingly modified 
[34l I35I . It has also been noted that the threshold for 
SF coincides with that for spectral broadening [3(|, an- 
other manifestation of the material nonlincarity in the 
form of a rapidly time varying self phase modulation 
(SPM). For powers just above P cr it the beam will con- 
tract spatially and the energy will move toward the back 
of the pulse and will eventually lead to pulse steepening 
and the formation of a shock wave at the trailing edge 
[T2I W\ |3S| . The intensities reached by the pulse give rise 
to efficient CE while pulse steepening may explain the 
observed higher conversion efficiencies for the anti-stokes 
(blue-shifted) components As the power is further 
increased (P ~ 1.5 — 2 • P cr it) the pulse will undergo tem- 
poral splitting [U due to the effect of normal GVD. 
Even higher powers (P > 3-P cr it) will lead to a pulse that 
evolves by pushing the energy toward the leading edge 
[371 l39l liol Elj . The peak intensities have been shown 
to be sufficient to excite a plasma with a large enough 
concentration as to even compensate the Kerr induced 
refractive index variation |36t |42j . The same plasma will 
eventually reduce the peak intensity through absorption, 
but for sufficiently high powers multiple pulse reforma- 
tion occurs with each new pulse rising from the trailing 
background power [3^, |4jj. Along side an evident CE 
it has been observed that a large part of the generated 
spectrum is actually emitted along the pulse propagation 
axis |43j . This has been explained by taking into account 
a possible Kerr-lensing effect induced by the pulse itself 



3 



[36( or possibly by the formation of an effective waveg- 
uide |44|. From this brief overview we can see that a 
large number of phenomena are involved, some occurring 
simultaneously while the presence of others depends on 
the particular input conditions. In the following sections 
we shall describe the experimental investigation method 
and results that shed some new light on the details of the 
filamentation process. 



III. EXPERIMENTAL SETUP 

We present data concerning filamentation in a 15 mm 
long, type II, lithium triborate (LBO) crystal. Such crys- 
tals are usually used for experiments involving the second 
order (x^) nonlinearity, however we rotated the crys- 
tal so that all second order processes are severely phase- 
mismatched (and thus negligible) whilst still maintain- 
ing the pump pulse propagation axis perpendicular to 
the input facet. The input laser pulse is delivered from 
a frequency doubled, 10 Hz Nd:glass mode-locked and 
regeneratively amplified system (Twinkle, Light Conver- 
sion) and has 200 fs duration and 527 nm central wave- 
length. Spatial filtering of the pulse before entering the 
sample is necessary to guarantee a uniform gaussian-like 
profile that will thus seed only one filament at a time (i.e. 
no spatial break-up into multiple filaments is observed). 
The beam has a FWHM of 5 mm diameter and is then 
focused onto the LBO crystal with a 50 cm focal length 
lens placed at 50.1 cm from the crystal entrance facet. 
The angular spectra of the filament have been detected 
by an imaging spectrometer (Oriel Instruments, 77250- 
M with a 1200 1/mm grating) placed after the crystal. 
The device reconstructs without distortion the entrance- 
slit at the output plane (the slit image is not on a curved 
surface, as would occur for conventional optics), with the 
different frequency components at different lateral posi- 
tions. By placing the entrance slit in the focal plane of 
a focusing lens, the angular frequency distribution can 
be detected in a single shot. The actual spectra are cap- 
tured by placing a CCD camera in the monochromator 
output imaging plane: the central regions require a very 
high dynamical range 16 bit CDD camera (Andor EEV 
40-11) in order to avoid saturation and loss of low-power 
details whereas the outer regions were characterized us- 
ing a higher spatial-resolution 8 bit camera (Pulnix TM-6 
CN). Note that only the angular distribution in the slit 
plane will be monitored, which is sufficient as long as the 
axial symmetry of the process is preserved. In acquir- 
ing the angular spectra it is very important to guarantee 
the possibility of single shot acquisition. In fact while 
the envelope spectral shape gives the information on the 
small-scale structure formed in the near field, it will be 
the fine interference-fringe structure in the spectra what 
reveals the extended feature in the space-time domain. 
For example, a single X-wave, a couple, a train and even 
a chaotic gas of X waves will all produce the same enve- 
lope angular spectrum, while the difference between the 
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FIG. 2: Two examples of pulse-filamentation far-field spectra 
taken for two different materials: (a) in 15 mm of water, 
P ~ 2 ■ P C rit and (b)in 15 mm of LBO, P ~ 1 • P cr it. (b) 
was obtained from three separate images each one spanning a 
different wavelength range (420-510 nm, 510-540 nm, 540-650 
nm) 

possible realizations will appear only in the modulated 
fine structure of the spectrum. These fine details will 
fluctuate slightly from shot to shot, owing to input-pulse 
energy, duration and diameter fluctuations. 



IV. EXPERIMENTAL RESULTS 

In figureQwe show the filament 1/e 2 diameter against 
input pulse energy. We see that as power increases the 
beam diameter contracts until it reaches a "stationary" 
value - we identify the lowest energy for which this 
occurs as the corresponding critical threshold energy, 
Eth- We therefore performed all our measurements in 
the energy range for which stable single filamentation 
occurs, i.e. from 2 /xJ to 5.6 //J. In order to illustrate that 
the phenomenum under inspection is a general feature 
of optical pulse propagation and is not particularly 
related to the specific material, we show in figure [3 two 
angular spectra, the first obtained in water (fig|2ja)) 
with E ~2-E t h and the second in LBO (Fig. |2») with 
E ~ 1 • E t h- Fig. Efb) was obtained from three separate 
images each one spanning a different wavelength range 
(420-510 nm, 510-540 nm, 540-650 nm) so as to minimize 
the effects of "blooming" from the central high-intensity 
region while still keeping significant detail in the lower 
intensity regions. Note also that very similar spectral 
features have also been reported elsewhere relative to 
filamentation of picosecond pulses in air [3(j . Although 
FigEIa) has fewer details due the large number (30) of 
shots over which the profile for water was integrated, the 
spectra have in common a definite "X-likc" pattern with 
a strong on-axis emission. We obtained similar figures 
also for fused silica. 

We shall now look more closely at the angular spectra, 
paying particular attention to the details. In figures 
|3Ja), (b) and (c) we show the angular spectra for input 
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FIG. 3: Stokes component of the measured far-field spectra, 
(a) Input pulse energy Ei„ = 2 fiJ: in the inset we show an- 
other spectrum taken for a slightly smaller Ei„ that highlights 
the modulation fringes. The total spectrum is showed in Fig 
El (b) Ein — 3 /iJ: in the inset we show the total recorded 
spectrum, (c) Ei n — 4.8 pJ. Note the different modulation 
patterns in the three measurements. 



energies E in equal to 2, 3 and 4.8 /xJ, respectively, 
thus mapping the evolution of the pulse structure for 
increasing powers. In all these figures we have focused 
our attention on the Stokes part of the spectrum: the 
complete recorded spectrum for E% n — 2 /iJ is shown 
in Fig. Efb) and that for Ei n = 3 /iJ in the inset of 
Fig. 02h). We first note the strong on-axis (i.e. for 
small transverse wave-vectors or angles) emission (AE) 
that extends well both into the blue Stokes and red 
anti-Stokes regions, characteristic of pulse filamentation 
and subsequent continuum generation. If considered in 
frequency rather than in wavelength, the anti-Stokes 
components have a larger extension (~ 0.9 fs _1 ) than 
the Stokes components (~ 0.7 fs _1 ), in agreement with 



literature 0, H^, ^5[. Together with the low-angle 
emission there is also a distinct "X" pattern, a signature 
of CE with a much more pronounced extension of the 
anti-Stokes component (~ 0.9 fs _1 to be compared with 
the ~ 0.2 fs" 1 of the Stokes part). All of the recorded 
spectra show a surprisingly regular pattern for small 8. 
This pattern is not so obvious in Fig. EJa) so we have 
included in the inset the Stokes spectrum for a different 
laser shot in which it is much clearer. The fringes are 
centered at 9 = with parabolic-like dependence on 
wavelength (similar features have also been observed 
in air with longer input pulse durations |30j|'). It is 
interesting to note how the sign of the fringe curvature 
inverts in passing from E; m — 2 fiJ to Bj„ = 3 fiJ and 
then shows a further sovra-modulation for Ei n = 4.8 /xJ 
(Fig. |3J|. We note that this behavior was found also in 
other materials (e.g. water) and under different focusing 
conditions with the only difference being the actual 
input energies at which the various modulation patterns 
are observed. 



V. DISCUSSION 

Despite the large amount of data available in litera- 
ture, to the best of our knowledge none of the numerical 
simulations shown to date display a combination of the 
main features we have measured. Namely, these are 
the distinct X-arms, the strong and largely extended 
axial emission, the periodic modulation of the axial 
emission and, finally, the inversion of this modulation 
pattern with increasing input power. Given this large 
difference, we try to give an explanation of the spectra 
by studying the properties of a linear superposition 
of stationary states, with which we approximate the 
istantaneous wave forms (into which the actual pulse 
shape may be decomposed) inside the material. The 
nature of the linear states will depend on the nature of 
the associated spectral shape so that CE will be related 
to an X-like wave [l4| whereas AE may be simply related 
to a gaussian-like spatio-temporal profile. We underline 
that the following discussion is not aimed at explaining 
beam filamentation but rather at trying to justify the 
experimental spectral features starting from results 
presented in literature. 

We start by analyzing the spectrum for Ei n = 2 /iJ. 
The input power is just above P cr it so that during 
propagation the pulse starts to collapse - in doing so a 
strong CE is initiated and finally the energy is moved 
to the trailing edge of the filament, as discussed above. 
Therefore we may describe the overall pulse as an 
X-wave (associated to the CE) and a trailing gaussian 
profile. The X wave is described following the recipe 
given elsewhere |32j with the input parameters taken 
as the material refractive index and dispersion relations 
and the total bandwidth of the actually measured CE 
(110 nm). The Gaussian profile is chosen so that it 
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has a 1/e 2 width equal to that measured (26 fim) and 
a temporal bandwidth corresponding to the measured 
AE bandwidth (230 nm). We note that all the results 
described below do not actually depend on the particu- 
lar function chosen to describe the gaussian-like profile 
and no substantial variations were observed using a 
super-gaussian or hyperbolic-secant profile with respect 
to a simple gaussian form. Furthermore we may expect 
a temporal delay between the two profiles and indeed 
it is the interference between these two that gives rise 
to the measured modulation patterns. However, we 
are still missing an ingredient. If we take the above 
described pulses, and Fourier-transform the sum of these 
so as to obtain the angular spectra we observe only 
straight, vertically aligned interference fringes. In order 
to explain the curvature of these we must also introduce 
a transverse spatial chirp (i.e. a phase front curvature) 
into the gaussian profile, ft is well known that in the 
presence of GVD and/or SPM the pulse may develop 
a strong temporal chirp. It has also been noted that 
in the simultaneous presence of SPM and SF, the pulse 
may also develop a strong frequency dependent mode 
size (FDMS) |4g, and a spatial chirp. In analogy with 
a linear temporal chirp |48} we may write the complex 
spatial amplitude as A(r) — exp[— (1 + ja r ) ■ (r/a) 2 ] 
were a is the beam width and a r is the spatial chirp 
parameter. Fig. EI a ) shows the (9 — A) distribution 
for an X wave and a spatially chirped gaussian pulse 
with a r = —2 and temporally delayed by r = —200 fs, 
thus approximating the power build-up at the trailing 
edge of the pulse expected for low input powers. The 
absolute value of a r was chosen so as to match the 
angular divergence of the AE and the sign so as to 
match the measured fringe curvature direction. Indeed, 
as can be seen the axial component shows a definite 
curved fringe pattern in close agreement with that 
shown in figure |3Ja). A negative spatial chirp (i.e. a 
defocusing phase-front curvature) could be due to a 
delayed-plasma induced defocusing. However we note 
that another picture is also possible and that is one in 
which the spatial chirp has opposite sign (a r = +2) 
and the gaussian pulse (shock wave) is in front of the 
now trailing X-wave. Although a leading shock front 
with a focusing wave-front is in fact compatible with a 
strong SF regime, we believe the picture of a trialing 
shock wave (at input powers near threshold) to be more 
acceptable in the frame of the present literature. 

Regarding the presence of a strong axial emission, the 
scientific community is still lacking a widely accepted 
theory although the most often used explanations involve 
the formation of an effective waveguide |44j induced by 
the balance of plasma defocusing and SF or just from 
the effect of SF alone j3||- However SF is just the spatial 
manifestation of SPM. The temporal analogue generates 
new frequencies and, in the spirit of this paper, the two 
should not be considered separately. Indeed, it is the 
spatio-temporal effect of SPM that explains CE, so it 
is not obvious why at a certain point of the pulse evo- 
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FIG. 4: Calculated far-field spectra of stationary linear 
states, (a) a gaussian pulse with spatial chirp a s = —2 and 
temporally delayed by +200 fs, with respect to an X-wave, (b) 
two equal-intensity spatially chirped gaussian pulses tempo- 
rally shifted with respect to an X-wave by -220 fs (q s = +2) 
and +220 fs (a s — —2) and (c) three equal-intensity spatially 
chirped gaussian pulses temporally shifted by -300 (a s = +2), 
-200 (a s = +2) and +300 fs (a s = -2) with an X-wave in 
t=0. All gaussian radial diameters are 25 pan at 1/e 2 . 



lution SPM should split into separate spatial (SF) and 
temporal (spectral broadening) phenomena. It could be 
possible that nonlinear losses (NLL) play an important 
role. Numerical simulations confirm that in the presence 
of NLL a high intensity pulse will develop a flat top (as 
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also reported in literature 47]) and, correspondingly, 
the spectral components generated via SPM will be 
strongly limited. This should occur in both the spatial 
and temporal coordinates. However there may exist 
mechanisms which break the spatio-temporal symmetry, 
such as Raman nonlinearity or plasma generation, which 
in turn could force the wave collapse in the sole temporal 
dimension, i.e. the formation of a temporal shock-wave 
with a spatial flat-top. In this case we would observe 
AE, an efficient broadening of the temporal spectrum 
with spatial wave-vectors aligned along the input pulse 
propagation axis. 

As for the spatial chirp, we note that the fringe pattern 
for higher powers is qualitatively different from that 
observed near P cr it- Indeed not only is the curvature 
inverted but the modulation is also much deeper and it 
is not concentrated near the central carrier frequency 
but rather becomes clearer at frequency shifts greater 
than those obtained by CE. If we take into account 
these features we realize that the interference is not due 
to an interplay between a gaussian and an X profile 
but rather between two (or more) gaussian-like pulses. 
Moreover, the chirp of the two pulses must have opposite 
sign and similar absolute value if such sharp and deep 
modulation patterns are to be explained. In Fig. Efb) 
we kept a r — — 2 for the rear pulse and put a r = +2 for 
the leading pulse. This condition gives us high contrast 
fringes with the same measured curvature shown in fig. 

Hb). 

The temporal delays of the two gaussian pulses in fig. 
0Jb) where chosen so as to match the experimental 
fringe frequency, so that r = —220 and r = 220 fs with 
respect to a central X-wave. As already discussed, we 
are to expect pulse splitting or dynamical replenishment 
for these higher input energies. The actual position 
of the X component does not induce any relevant 
changes in the AE and we are not able to precisely 
determine the temporal location of the CE. This may 
be explained by noting that the total power contained 
in the experimental spectrum related to CE part (after 
the monocromator input slit) is always at least a factor 
ten smaller than that due to AE, so that interference 
between the two contributions is rather weak (see for 
example the very low contrast interference fringes in figs. 
E^a) andQJa)). The total delay between the leading and 
trailing pulses (440 fs) is surprisingly high if compared 
to the 200 fs input pulse duration. Obviously GVD 
is playing a major role and is indeed dominating the 
temporal profile evolution. 

Figure 0Jc) was obtained with three such gaussian 
pulses with temporal delays r = —400, —240, +400 and 
a = +2, +2 and —2 respectively. Once again we were 
not able to precisely position the X-wave as negligible 
variations were observed. The agreement with the exper- 
imental data (Fig. EUc)) is very good and the model is 
in agreement with other reported findings relative to the 
formation of multiple peak re-formation for high enough 
input energies [4(J, H|| ■ Furthermore, the temporal delay 



between the pulses has increased to 800 fs, i.e. four times 
the input pulse duration and indicates that for higher 
powers the interplay between SPM, SC generation and 
GVD is further enhanced, as expected. This delay is 
extremely large if compared to the input pulse duration 
however a further indication that this result is correct is 
given by the multiple-shot autocorrelation trace shown 
in fig. [5] obtained by imaging the output crystal facet 
onto the autocorrelation nonlinear crystal. The five 
peaks in the trace show that with a similar average 
pulse input energy we have multiple pulse splitting with 
the formation of three daughter pulses and a temporal 
separation of 600 fs between the leading and trailing 
pulses. 

So far we have purposely neglected temporal chirp 
(at) which is expected to be at least as important as 
the spatial chirp. However AE does not seem to be 
sensitive to this parameter. Indeed, variations of a t 
did not produce a significant change in the numerical 
AE fringe pattern but rather only a local reduction of 
contrast, the position of which depends on the sign and 
value of at in the various linear components with which 
the filament was modeled. We note though, that the 
experimental spectra also show a modulation in the CE 
in the form of a multiple "X-arm" splitting. In fact, by 
slowly increasing the pulse input energy from below to 
above threshold we never experimentally observed the 
formation of a single X pattern - CE emission in our 
experimental setup always seems to appear in the form 
of multiple X-like patterns. This sheds some light on the 
nature of the CE in the sense that it may be explained 
by super-imposing multiple X pulses so that the number 
of arms corresponds directly to the number of X pulses 
and each X must be spectrally shifted with respect to 
the others. This is physically feasible in the presence 
of a single input pulse that undergoes sever temporal 
chirping and then splits into multiple "X-pulses" (i.e. 
daughter pulses that gives rise to CE) so that each of 
these is centered on a different carrier frequency. The 
values for the spectral shift observed in our spectra vary 
from 0.05 to 0.2 fs -1 with the higher values observed at 
higher input energies. 



VI. CONCLUSIONS 

In conclusion we have shown the possibility of recov- 
ering detailed information regarding space-time coupled 
phenomena from angular-spectral characterization. The 
nature of the numerical model used to interpret the ex- 
perimental measurements, based on a linear combination 
of stationary states rather than on a full nonlinear evo- 
lution simulation, allows only a qualitative guess at the 
parameters involved. Nevertheless the simplicity of this 
approach is very appealing and indeed, by using a com- 
bination of linear states it is possible to understand if the 
pulse has changed its spatio-temporal energy distribution 
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FIG. 5: Autocorrelation trace of of the ouput pulse filament 
clearly showing that the pulse has split into three sub-pulses 
with a delay of 600 fs between the most external ones. E in ~ 4 

without actually experiencing splitting or if it has gone 
through a single or even multiple pulse reformation. The 
details of the spectra also reveal that the spatio-temporal 
coupling manifests itself in a strong spatial chirp of the 
filament. If the far- field measurements were to be com- 
pleted with a complete phase characterization we would 
of course be able to reconstruct the full space-time pro- 
file of the pulse through a Fourier transform - the lack of 
this information is reflected in the difficulty in retrieving 



more precise information for example of the exact tem- 
poral delay of the CE with respect to the AE sources. 
Furthermore, the spatial chirp parameter may be deter- 
mined precisely with the combination of an imaging spec- 
trometer and a spatial shearing interferometer [5fJ] and 
shall be considered in future measurements. We further 
underline that whereas the "X" part of the spectra are 
numerically well reproduced in the frame of our simple 
linear model in both the Stokes and anti-Stokes regions, 
this is not true for the AE spectral components. The anti- 
Stokes AE is rather confused and the (not always visible) 
fringe pattern seems to show a fast modulation frequency 
that is not compatible with the parameters that fit the 
Stokes region. This could indicate that the shock waves 
that generate the blue and red shifts in the spectrum are 
spatially distinct and have different temporal delays. 
Notwithstanding these shortcomings, the aim of this 
work is to underline the importance of angular spectral 
measurements that, together with other methods such as 
FROG or 3D-mapping, may give a complete and exhaus- 
tive characterization of nonlinear optical-wave collapse 
phenomena. 

** Present address, Institut de Ciencies Fotoniques, 
c/Jordi Girona 29, NEXUS II, E-08034 Barcelona, Spain. 
*** Present address, Laboratoire Kastler Brossel, Uni- 
versite P. ct M. Curie, Case74, 4, place Jussieu, F75252 
Paris, France 



[1] G. A. Askarjan, "Interaction of laser radiation with vi- 
brating surfaces," JETP, 42, 1672 (1962). 

[2] W. Talanov, Izv. Vysshikh Uchebn. Zavedenii, Ra- 
diofizika 7 (1964). 

[3] R.Y.Chiao, E. Garmire, and C.H.Townes, "Self-trapping 
of optical beams," Phys. Rev. Lett. 13, 479-482 (1964). 

[4] M. Hercher, J. Opt. Soc. Am. 54, 563 (1964). 

[5] A. Braun, G. Korn, X. Liu, D. Du, J. Squier, and 
G.Mourou, "Self-channeling of high-paek-power fem- 
tosecond laser pulses in air," Opt. Lett. 20, 73-75 (1995). 

[6] K. Wilson and V. Yakovlev, "Ultrafast rainbow: tunable 
ultrashort pulses from a solid-state kilohertz system," J. 
Opt. Soc. B 14, 444-448 (1997). 

[7] R. Alfano, The Supercontinuum Laser Source (Springer- 
Verlag, New York, 1989). 

[8] P. Rairoux, H. Schillinger, S. Niedermeier, M. Rodriguez, 
F. Ronneberger, R. Sauerbrey, B. Stein, D. Waite, C. W. 
nd H. Wille, L. Woste, and C.Ziener, "Remote sensing 
of the atmosphere using ultrashort laser pulses," Appl. 
Phys. B 71, 573-580 (2000). 

[9] P. Nozieres and D. Pines, The Theory of Quantum Liq- 
uids, vol. II (Addison- Wesley, Redwood City, 1990). 
[10] D.Strickland and P. Corkum, "Resistance of short pulses 

to self-focusing," J.Opt.Soc.Am.B 11, 492-497 (1994). 
[11] G. Luther, A. Newell, J. Moloney, and E. Wright, "Short 
pulse conical emission and spectral broadening in nor- 
mally dispersive media," Opt. Lett. 19, 789-791 (1994). 
[12] A. Gaeta, "Catastrophic collapse of ultrashort pulses," 



Phys. Rev. Lett. 84, 3582-3585 (2000). 

[13] J. Ranka and A. Gaeta, "Breakdown of the slowly varying 
envelope approximation in the self-focusing of ultra-short 
pulses," Opt. Lett. 23, 534-536 (1998). 

[14] S. Trillo, C. Conti, P. D. Trapani, O. Jedrkiewicz, 
J. Trull, G. Valiulis, and G. Bellanca, "Coloured coni- 
cal emission via second-harmonic generation," Opt. Lett. 
27, 1451-1453 (2002). 

[15] G. Valiulis, J. Kilius, O. Jedrkiewicz, A. Bramati, S. Mi- 
nardi, C. Conti, S. Trillo, A. Piskarskas, and P. Di Tra- 
pani, "Space-time nonlinear compression and three- 
dimensional complex trapping in normal dispersion,"" in 
OSA Trends in Optics and Photonics (TOPS) 57, Quan- 
tum Electronics and Laser Science Conference (QELS 
2001), Technical Digest, Post conference Edition (Optical 
Society of America, Washington DC, 2001), pp. QPD10 
-1-2. 

[16] P. Di Trapani, G. Valiulis, A. Piskarskas, O. Jedrkiewicz, 
J. Trull, C. Conti, and S. Trillo, "Spontaneously gen- 
erated X-shaped light bullets," Phys. Rev. Lett. 91, 
093 904-1 (2003). 

[17] C. Conti, S. Trillo, P. Di Trapani, G. Valiulis, 
A. Piskarskas, O. Jedrkiewicz, and J. Trull, "Nonlinear 
electromagnetic X waves," Phys. Rev. Lett. 90, 170 406- 
1 (2003). 

[18] O. Jedrkiewicz, J. Trull, G. Valiulis, A. Piskarskas, 
C. Conti, S. Trillo, and P. Di Trapani, "Nonlinear X 
wavesin second harmonic generation: esperimental re- 



8 



suits," Phys. Rev. E 68, 026 610-1 (2003). 

[19] D.Salerno, O. Jedrkiewicz, P. Di Trapani, J. Trull, and [35 
G. Valiulis, "Impact of dimensionality on noise-seeded 
modulational instability," in Nonlinear guided waves and 
their applications (2004), pp. WA-7. [36 

[20] S. Minardi, J. Yu, G. Blasi, A. Varanavicius, G. Valiulis, 
A. Berzanskis, A. Piskarskas, and P. Di Trapani, "Red 
solitons: evidence of spatiotemporal instability in x' 2 ' [37 
spatial soliton dynamics," Phys. Rev. Lett. 91, 12 390-1 
(2003). 

[21] M. Kolesik, E. Wright, and J. Moloney, "Dynamic Non- [38 
linear X-waves for Femtosecond Pulse Propagation in 
Water," Arxiv:Phys/0311021vl (2003). 

[22] A. Dubietis, E. Gaizauskas, G. Tamosauskas, and 

P. Di Trapani, "Light filaments without self-channeling," [39 
accepted for publication in Phys. Rev. Lett. (2004). 

[23] S. Diddams, H. Eaton, A. Zozulya, and T. Clement, 

"Amplitude and phase measurements of femtosecond [40 
pulse splitting in nonlinear dispersive media," Opt. Lett. 
23, 379-381 (1998). 

[24] A. Bernstein, J. Diels, T. Luk, T. Nelson, A. McPherson, [41 
and S. Cameron, "Time resolved measurements of self 
focusing pulses in air," Opt. Lett. 28, 2354-2356 (2003). 

[25] L. Gallmann, G. Steinmeyer, D. Sutter, T. Rupp, C. Ia- 

conis, I. Walmsley, and U. Keller, "Spatially resolved am- [42 
plitude and phase characterization of femtosecond optical 
pulses," Opt. Lett. 26, 96-98 (2001). 

[26] C. Dorrer, E. Kosik, and I. Walmsley, "Direct space-time 

characterization of the electric fields of ultrashort optical [43 
pulses," Opt. Lett. 27, 548-550 (2002). 

[27] M. Potenza, S. Minardi, J. Trull, G. Blasi, D. Salerno, 
A. Varanavicius, A. Piskarskas, and P. Di Trapani, 
"Three dimensional imaging of short pulses," Opt. Com- [44 
mun. 229, 381-390 (2004). 

[28] J. Trull, O. Jedrkiewicz, P. D. Trapani, A. Matijosius, 
A. Varanavicius, G. Valiulis, R. Danielius, E. Kucinskas, 
and A. Piskarskas, "Spatiotemporal three-dimensional [45 
mapping of nonlinear X-waves," Phys. Rev. E 69, 
026 607-1 (2004). 

[29] A. Matijosius, J. Trull, P. Di Trapani, A. Dubietis, 

R. Piskarskas, A. Varanavicius, and A. Piskarskas, "Non- [46 
linear space-time dynamics of ultrashort wave packets in 
water," Opt. Lett. 29, 1123-1125 (2004). 

[30] D. Mikalauskas, A. Dubietis, and R. Danielius, "Observa- [47 
tion of light filaments induced in air by visible picosecond 
laser pulses," Appl. Phys. B 75, 899-902 (2002). 

[31] A. Dubietis, G. Tamosauskas, I. Diomin, and [ 
A. Varanavicius, "Self-guided propagation of femtosec- 
ond light pulses in water," Opt. Lett. 28, 1269-1271 
(2003). [49 

[32] M. Porras and P. D. Trapani, "Localized and station- 
ary light wave modes in dispersive media," Los Alamos 
Arxive physics/0309084 (2003). 

[33] J. Marburger, "Self-focusing: theory," Progr. Quant, [50 
Electron. 4, 35-110 (1975). 

[34] P. Chernev and V. Petrov, "Self-focusing of light pulses 
in th presence of normal group- velocity dispersion," Opt. 



Lett. 17, 172-174 (1992). 

G. Luther, J. Moloney, A. Newell, and E. Wright, "Self- 
focusing threshold in normally dispersive media," Opt. 
Lett. 19, 862-864 (1994). 

A. Brodeur and S. Chin, "Ultrafast white light contin- 
uum generation and self-focusing in transarent condensed 
media," J.Opt.Soc.Am.B 16, 637-650 (1999). 
A. Zozulya and S. A. Diddams, "Dynamics of self-focused 
femtosecond laser pulses in the near and far fields," Opt. 
Expr. 4, 336-343 (1999). 

J. Rothenberg, "Space-time focusing: breakdown of 
the slowly varying envelope approximation in the self- 
focusing of femtosecond pulses," Opt. Lett. 17, 1340- 
1342 (1992). 

M. Mlejnek, E. Wright, and J. Moloney, "Dynamic spa- 
tial replenishment of femtosecond pulses propagating in 
air," Opt. Lett. 23, 382-384 (1998). 

M. Mlejnek, E. Wright, and J. Moloney, "Power depen- 
dence of dynamic spatial replenishment of femtosecond 
pulses propagating in air," Opt. Expr. 4, 223-228 (1999). 
O. Kosareva, V. Kandidov, A. Brodeur, C. Chen, and 
S. Chin, "Conical emission from laser-plasma interactions 
in the filamentation of powerful ultrashort laser pulses in 
air," Opt. Lett. 22, 1332-1334 (1997). 
S. Tzortzakis, L. Sudrie, M. Franco, B.Prade, A. Mysy- 
rowicz, A. Couairon, and L. Berge, "Self-guided propa- 
gation of ultrashort IR laser pulses in fused silica," Phys. 
Rev. Lett. 87, 3902-1 (2001). 

W. Smith, P. Liu, and N. Bloembergen, "Superbroaden- 
ing in H2O and D2O by self focused picosecond pulses 
from a YAlG:Nd laser," Phys. Rev. A 15, 2396-2403 
(1977). 

E. Nibbering, P. Curley, G. Grillon, B. Prade, M. Franco, 

F. Salin, and A.Mysyrowicz, "Conical emission from self 
guided femtosecond pulses in air," Opt. Lett. 21, 62-64 
(1996). 

W. Liu, O. Kosareva, I. Golubtsov, A. Iwasaki, 
A. Becker, V. Kandidov, and S. Chin, "Femtosecond laser 
pulse filamentation versus optical breakdown in H2O" 
Appl.Phys. B 76, 215-229 (2003). 

S. Cundiff, W. Knox, E. Ippen, and H. Haus, "Frequency 
dependent mode size in broadband Kerr-lens mode lock- 
ing," Opt. Lett. 21, 662-664 (1996). 
S. Polyakov, F. Yoshino, G. Stegeman, "Interplay be- 
tween self-focusing and high-order multiphoton absorp- 
tion," J. Opt. Soc. Am. B 18, 1891-1895 (2001). 
S. Akhmanov, V. Vysloukh, and A. Chirkin, Optics 
of Femtosecond Laser Pulses (American Institute of 
Physics, New York, 1992). 

W. Liu, S. Chin, O. Kosareva, I. S. Golubtsov, and 
V. Kandidov, "Multiple refocusing of a femtosecond laser 
pulse in a dispersive liquid (methanol)," Opt. Commun. 
225, 193-209 (2003). 

C. Dorrer and I. Walmsley, "Simple linear tecnique for 
the measurement of space-time coupling in ultrashort op- 
tical pulses," Opt. Lett. 27, 1947-1949 (2002). 



